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A key point in the spin foam approach to quantum gravity is the implementation of simpUcity 
constraints in the partition functions of the models. Here, we discuss the imposition of these 
constraints in a phase space setting corresponding to simplicial geometries. On the one hand, this 
could serve as a starting point for a derivation of spin foam models by canonical quantisation. On 
the other, it elucidates the interpretation of the boundary Hilbert space that arises in spin foam 
models. 

More precisely, we discuss different versions of the simplicity constraints, namely gauge-variant 
and gauge-invariant versions. In the gauge-variant version, the primary and secondary simplicity 
constraints take a similar form to the reality conditions known already in the context of (complex) 
Ashtekar variables. Subsequently, we describe the effect of these primary and secondary simplicity 
constraints on gauge-invariant variables. This allows us to illustrate their equivalence to the so- 
called diagonal, cross and edge simplicity constraints, which are the gauge-invariant versions of the 
simplicity constraints. In particular, we clarify how the so-called gluing conditions arise from the 
secondary simplicity constraints. Finally, we discuss the significance of degenerate configurations, 
and the ramifications of our work in a broader setting. 



I. INTRODUCTION 

In the quest to find a theory of quantum gravity, a class of theories arise, whose respective Hilbert spaces 
are labelled by graphs. This class includes loop quantum gravity spin foams quantum Regge calculus 
3, causal dynamical triangulations Q and so on. In loop quantum gravity, which is a canonical quantisation 
of continuum general relativity, these graphs (or spin networks as they are known) arise naturally and provide 
a basis for the space of states prior to the imposition of the Hamiltonian constraint. On the contrary, in spin 
foams and quantum Regge calculus, the graphs take on a greater significance, implementing a discretisation 
of the continuum data. In these quantisation methods, one views the abstract graph as dual to a simplicial 
manifold, upon which the gravitational theory is expressed.^ Thus, it acts as a uv-regulator which naturally 
cuts out certain divergences from the quantisation. There is an obvious trade-off here. Altering the fundamental 
component of space-time, the manifold itself, in a theory meant to describe the quantum dynamics of that very 
structure means that one must tackle anew some subtle technical issues. 

Already at the classical level, one needs to formulate a consistent lattice dynamics, and furthermore, one should 
understand the significance and extent to which the lattice structure affects a realisation of the symmetries of 
the theory 0, Q . Even before dealing with this issue, upon a simplicial manifold, the ease with which one can 
formulate an action principle and endow the topological structure with gravitational data depends heavily upon 
how one parameterises the theory. A natural way to encode a simplicial metric geometry is to use edge-lengths. 
Indeed, this is the underlying structure of Regge's original proposal. These are not the variables one uses, 
however, in the framework upon which modern spin foam quantisation rests. 

Spin foam quantisation is implemented on a class of theories that emerge as the lattice counterparts of non- 
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metric formulations of general relativity, in particular the Plebanski formulation 9]. In the continuum, these 
actions are functionals of a bi-vector, a spin connection and Lagrange multipliers that ensure the solution space 
includes that of the metric theory.^ These actions take the form of constrained topological field theories, where 
the unconstrained action describes BF theory. The advantage is that one circumvents initial stumbling blocks 
which plague other formulations. The path integral quantisation of BF theory on a simplicial manifold is well 
understood. The measure on the space of paths is well-defined and anomaly-free while the weight for each path 
is given succinctly by a lattice form of the action. Since any quantum amplitude is a topological invariant, one 
may view it as a well-defined prescription for the amplitude on any simplicial or continuum manifold of the same 
topology. The aim is to devise a similar strategy for the constrained theory. First one must devise a system of 
classical constraints for a simplicial manifold which restrict the topological theory to the gravitational theory. 
Then, one must impose those constraints consistently at the quantum level. This will alter both the measure 
and the weights of the original topological amplitude. 

The process is non-trivial for a multitude of reasons and constitutes the key point in the construction of spin 
foam models. Different methods to impose the constraints in the discrete setting have been developed and in 
fact, this partially constitutes the differences among these existing models [Tol - [l5| . In a canonical analysis of the 
continuum Plebanski action, both primary and secondary simplicity constraints arise. The secondary set ensure 
that the primary set are preserved under time-evolution. Most of the aforementioned methods implement only 
a discrete version of the primary simplicity constraints explicitly in the path integral amplitudes. In contrast, 
[l6| argues that one should impose both the primary and secondary simplicity constraints, as this would ensure 
the correct correlation functions for all observables. This proposal has not been brought to completion yet, as 
a suitable discretisation for the secondary simplicity constraints was missing [l^. Additionally, the secondary 
simplicity constraints involve the connection variables, which are far more difficult to constrain within the spin 
foam path integral than the bi-vectors. 

In this work, we shall discuss a phase space corresponding to a discretised Plebanski theory, in particular we 
shall develop discrete versions of the primary and secondary simplicity constraints. One motivation is that the 
construction of a phase space path integral for the discrete theory, as suggested in [l^, clearly needs such an 
analysis. In particular, for the phase space path integral measure, one needs to find the symplectic structure 
on the reduced phase space. 

Another reason arises upon finding that current spin foam models encode a boundary Hilbert space [l^ 
that is related to the Hilbert space of loop quantum gravity restricted to a particular graph (i.e. a discrete 
version of loop quantum gravity). One should examine, however, the discussion in Jl6J expressing an alternative 
viewpoint. Clearly, this boundary Hilbert space should be related to a quantisation of a phase space reduced by 
the simplicity constraints. As was already pointed out in (2]| . once one imposes all the simplicity constraints, the 
reduced phase space is much smaller than the phase space corresponding to a discretised version of loop quantum 
gravity. It is important to understand this discrepancy in more detail and we shall comment on this issue in 
the course of this work. Also, if one were to derive a spin foam model from a canonical quantisation, following 
a standard implementation of the (second class) simplicity constraints, one would improve the reliability of the 
spin foam approach. The standard way to enforce second class constraints is to solve them classically. Therefore, 
we need a detailed knowledge of the reduced phase space. A canonical realisation of the discrete theory in the 
Plebanski formulation and its subsquent quantisation would also enlighten the connections between spin foam 
models and loop quantum gravity. 

A canonical analysis of Plebanski theory in the continuum is provided by ,1^] . The analysis of the discrete 
case is usually much more subtle. There is no guarantee that the processes of Legendre transformation and 
discretisation commute. In particular, discretisation breaks diffeomorphism symmetry 0,0- The Hamiltonian 
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and diffeomorphism constraints, which determine the canonical dynamics in the continuum, follow from the 
invariance of the action under diffeomorphisms. Therefore, the breaking of diffeomorphism symmetry in the 
discrete setting alters the concepts of the constraints; for more details, see 0, Q- This is one of the reasons why 
so far there is not a consistent (i.e. first class) set of discretised Hamiltonian and diffeomorphism constraints 
for 4d gravity. Naturally, this affects the computation of the secondary simplicity constraints, which in the 
continuum are determined by the time evolution of the primary constraints. More precisely, they are defined 
as the Poisson brackets between these primary constraints and the Hamilltonian constraints. Thus, we cannot 
derive the secondary constraints along conventional lines. 

To sidestep such difficulties, we proceed along a different route that was also followed in [2l| - |23j . The main 
task is to identify a set of constraints on the phase space data (bi- vectors and holonomies of the spin connections) 
that spares only geometric configurations. In other words, the bi-vectors should be constructible from tetrads, 
while the spatial part of the spin connection should be the Levi-Civita connection. Such constraints were 
proposed in I21I-I23I . However, the works [l^,!!^! did not analyse the reduced phase space, which has now been 
performed in (2ll|. There, it was found that reducing the phase space does not correspond to just replacing the 
gauge group S0(4) with S0(3), which is what one might first expect if one considers the primary simplicity 
constraints only. Rather, if one takes all the constraints into account, one obtains a much smaller phase space. ^ 
In this work, we shall unravel the underlying mechanism. The basic idea is the following. 

The (canonical) Plebanski formulation starts from so(4)-valued bi-vectors Xij^ and S0(4)-valued holonomies 
Mij^^ , both of which are associated to the triangles {ij} of the 3d triangulation. These variables are subject 
to the diagonal and cross simplicity constraints - which correspond to the primary simplicity constraints - as 
well as the so-called edge simplicity constraints [2T'| - which correspond to the secondary simplicity constraints. 
Additionally, there are the (first class) Gaufi constraints arising from the SO (4) gauge symmetry. 

The diagonal simplicity and cross simplicity constraints are well known from the covariant discretisation of 
the Plebanski formulation [25!] . The geometrical meaning of these constraints transpires from the decomposition 
of so(4) into self-dual and anti-self-dual sectors: so(4) = su(2)+ ©su(2)_. The diagonal and cross simplicity 
constraints impose that areas and 3d-dihedral angles constructed in the self-dual sector coincide with those con- 
structed in the anti-self-dual sector. As these quantities allow one to reconstruct the 3d geometry of tetrahedra, 
one can summarise these constraints as requiring that 'left-handed geometry equals right-handed geometry' 
for each tetrahedron, as in the title of [2^ . One could therefore assume that the imposition of the full set of 
constraints determines the variables in one of the two chiral sectors as functions on the other sector. 

Indeed the secondary (or edge) simplicity constraints may also be cast in a form which imposes the equality 
of left-handed and right-handed quantities. But to be able to incorporate all these equalities, conditions have 
to be satisfied involving quantities of one and the same chiral sector. Importantly, they actu ally incorporate 
conditions among quantities within one chiral sector. These are the gluing constraints [2l|, [23I that ensure 
that the tetrahedra can be consistently glued along their common triangles. The gluing constraints are implied 
by the secondary simplicity constraints for the following reason. The secondary simplicity constraints require 



^ A way out is to consider the canonical analysis of the discrete action as is done in 0, 01 • The act of discretisation breaks 
the symmetries of the action and results in the replacement of constraints by pseudo-constraints. These pseudo-constraints are 
equations of motion, in which the canonical data of two consecutive time steps are only very weakly coupled to each other 
(Remember that true constraints are equations of motion, in which the canonical data of consecutive time steps are not coupled 
to each other at all. True constraints concern only the data of a single time step.) For discrete gravity one can typically only 
obtain an implicit version of these pseudo-constraints Q]. Moreover, such an analysis has so far not been performed for the 
Plebanski formulation. 

* This might be related to the difficulties encountered in [23|| in finding a real symplectic structure starting with the complex phase 
space of the self-dual Ashtekar formulation [23 | . Indeed, we shall find a form of the simplicity constraints resembling the reality 
conditons for the complex- valued Ashtekar variables. 
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that the spatial part of the S0(4) connection is of Levi-Civita type. Hence, it is determined by the spatial 
geometry, i.e. the bi-vectors. (In complex Ashtekar variables, which are related to the variables occurring in 
the self-dual Plebanski formulation, these constraints are the reality conditions stating that the real part of 
the Ashtekar connection is the Levi-Civita connection.) However, the spatial Levi-Civita connection, i.e. the 
rotations between the triads associated to neighboring tetrahedra, can be only determined if one is able to 
consistently glue tetrahedra together. Therefore, the gluing conditions must hold. 

Furthermore, the time components of the SO (4) connection encode the extrinsic curvature, which within a 
simplicial context is assimilated into the 4d-dihedral angles between neighboring tetrahedra. As we shall notice 
in the main text, the definition of the 4d-dihedral angles, using bi-vectors and holonomies, is somewhat abstruse. 
The gluing conditions ensure that the 4d-dihedral angles are unambiguously determined. The remaining part 
of the secondary simplicity constraints ensure that the 4d-dihedral angles constructed in the left-handed sector 
match those constructed in the right-handed sector, that is, that the self-dual and anti-self-dual geometries 
coincide. 

In Section [TTl we shall provide the essential foundations of a canonical formulation of discretised Plebanski 
theory. We shall proceed, in Section IIIIl by invoking time-gauge and deriving a set of constraints that ensure 
geometricity. In that setting, there is a nice separation of duties. One can clearly see that the primary simplicity 
constraints restrict the bi-vectors while the secondary constraints give conditions on the connection variables. 
To be able to satisfy all these secondary conditions further relations involving the bi-vectors must hold.^ 

Following that, in Section |IV[ we shall introduce gauge-invariant variables, namely the areas, 2d-, 3d-, and 
4d-dihedral angles and analyse the implications of the constraints, which were derived in the previous section, for 
these gauge-invariant variables. Moreover, we shall show that the diagonal, cross and edge simplicity constraints, 
as formulated in [2l|, , follow from the gauge- variant form of the constraints developed in Section IIIIl This 
proves the sufficiency of the gauge-variant primary and secondary constraints, since [21] contains a sufficiency 
proof of the diagonal, cross and edge simplicity constraints (provided the configurations are non-degenerate). 
In section |V] we shall point out the importance of non-degeneracy conditions. We shall see that without these 
conditions certain formulations of the simplicity constraints allow for a S0{3)-BF like sector of non-geometric 
solutions, whose dimension is larger than the one describing geometric configurations. We end with a discussion 
in Section IVll 



II. PASSING FROM THE CONTINUUM TO THE SIMPLICIAL 



In this section, we shall summarise the necessary fundamental facets of the discrete canonical formulation of 
4d BF and M Plebanski theory [H [ll]. 
The action in the continuum is: 



J M L 7 ^ 



(1) 



where X is the bi-vector, w is the spin connection, F is its curvature, and 7 is the Barbero-Immirzi parameter. 
The index A = aais anti-symmetric in a and a, with a, a £ {0, . . . , 3}, while A stands for the wedge product on 



In the discrete setting, the distinction between degrees of freedom encoded in the bi-vectors and degrees of freedom encoded in 
the connection or holonomies is not as clear-cut as in the continuum. Firstly, the definition of the discrete bi-vectors as continuum 
quantities integrated over the triangles in the triangulation involves the connection, see for instance [23l . Furthermore we shall 
adopt here a doubling of phase space variables, that is to every triangle we associate two bi-vectors and two holonomies. Between 
these pairs so-called matching conditions have to hold, that also mix the holonomies and the bi-vectors. See also the discussion 
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forms. The first term in the action describes a topological theory known as BF theory. The Lagrange multiplier 
^AB jj^^poses the simplicity constraints: 

X^AX^ = Ve^^, (2) 

where V ~ j^e^^ A . The non-degenerate class of solutions to these constraints falls into two sectors. 
More precisely, if V ie^^.p X^^ X^^ ^ 0, then 

4 r ±(eAe)^, 
^ = , . : A (3) 



± * (e A e) 

where e° is the co-frame 1-form field. These solution sets are referred to as topological and gravitational, 
respectively. For 7 = oo, substituting the topological solutions back into the action ((TJ produces an action that 
may be written as a total divergence, while if one substitutes for the gravitational solution set, one arrives at the 
Palatini action for general relativity. Of course for finite Immirzi parameter, such distinctions are not strictly 
valid as all the solutions give a mix of topological and gravitational actions, but one keeps the nomenclature in 
any case. 

In the continuum canonical analysis, one finds that there are primary simplicity constraints which are func- 
tionals of the bi- vectors only. These are not preserved under time-evolution unless certain secondary simplicity 
constraints are imposed. These secondary constraints are functions of both the bi-vectors and the connection 
variables. For those interested, we have summarised some properties of this canonical analysis in Appendix [Cl 
in particular, the form of the primary and secondary simplicity constraints. 

To obtain a discrete canonical formalism we intend to pursue a hybrid continuum-simplicial theory here, which 
allows for a straightforward construction of a symplectic form on phase space [13]. To do so, we keep time as a 
smooth parameter and assume that A4 can be folitated with leaves St, that is, = E x R. We triangulate this 
3d-hypersurface with flat tetrahedra: E — > A. The continuum fields are replaced by distributional analogues 
and integrated over suitable sub-simplices. 

Let us consider the BF theory alone for a moment. The canonical form of the action on A x R is: 



(4) 



where i,j label tetrahedra. Thus, the subscript {ij} picks out a triangle (since they are each shared by 
two tetrahedra). Xij and Mij are 
respectively. The tensor C"^^ is 1 
enforcing the following constraints: 



two tetrahedra). Xij and Mij are the bi- vector field and the parallel transport matrix between tetrahedra, 
respectively. The tensor C"^^ is the structure constant of so(4), while {1} is the set of Lagrange multipliers 



Irr A Y ^ M Y B 

^ij . — -^ji ' ji ij f 

:= E,^./, 



(5) 



where {1^1^+1} are the triangles sharing the edge. To each fiat tetrahedron, we associate a reference frame. 
Therefore, a priori, there are two bi-vectors associated to a triangle: Xij and Xji. The first constraint states 
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that they are related by an SO (4) transformation (that which parallel transports between the reference frames 
of the tetrahedra {i} and {j}). The second constraint states that the S0(4) transformation from tetrahedron 
{i} to tetrahedron {j} is the inverse of that from {j} to {i}. Thus, the first two constraints arise in the discrete 
formulation where we assign two pairs of canonical variables per triangle Xij , Mij and Xji , Mji . They ensure 
that these two pairs encode the same data (by expressing Xji, Mji as functions of Xij,Mij). Thus, we shall 
refer to them collectively as 'matching constraints'. The final two constraints are the closure (or GauB) and 
flatness constraints of lattice BF theory, respectively. 

We shall dispense with the flatness constraint since we shall not deal with the dynamical properties of gravity, 
let alone BF theory. We shall keep the closure constraint, since it is the lattice analogue of the Gaufi constraint 
that occurs in the canonical analysis of the Plebanski action (jC4l) . It generates S0(4) gauge transformations. 

It is convenient to separate the degrees of freedom into self-dual and anti-self-dual components using the 
Hodge-star operator on so (4). More details can be found in Appendix El The first term in the action tells us 
about the canonical commutation relations in the initial phase space, Vai upon which the constraints must act. 
The canonically conjugate variables are: {Mij±, Xij±), with commutations relations: 

7 ± 1 

7 ± 1 (6) 

{Xij±^,Xij±^} = ^ _^ ^(jABC Xij±'~' , 

{M,,±^^,M,,±'^^} = 0. 

The third relation is somewhat surprising, since it does not follow directly from the action, but it is necessary 
so that the commutators satisfy the Jacobi identity [1^ [s^l • The splitting into self-dual and anti-self-dual is 
respected by the symplectic structure, i.e. for any two combinations 1+ and of variables from the self-dual 
and anti-self-dual sectors respectively, we have {! + , Z_} = 0. 



Our aim is then to formulate a minimal set of constraints on the phase space Va , which remove non-geometrical 
canonical configurations. There are several ways to test whether such constraints do the job. But we shall call 
the constraints satisfactory if they allow us to reconstruct uniquely a vector for each edge of the 3d-hypersurface 
A. In other words, one can construct a co-frame for A from the Xij and Mij. This is exactly what the simplicity 
constraints accomplish in the continuum. Moreover, we might take some inspiration from the continuum theory 
and formulate the constraints as functions with a quadratic X dependence. Indeed, we adopted this form of 
the constraints in 1211 1: 





:- Xij+ 


■ Xij+ 


- Xij^ 


■ Xij-, 


&ljk 


:= Xij+ 


■ Xik+ 


- Xij^ 


■ Xik- , 


&ijkl 


:— Xik+ 




^Xjl+) 


— Xik- ■ {Mij-Xji- 



For the last constraints the triangles {ij}, {ifc}, {jl} must share an edge. These constraints are known as the 
diagonal, cross, and edge simplicity constraints. Although, we do not have a lattice Hamiltonian, we refer to the 
first two as primary, since they are independent of Mij , while we refer to the final one as a secondary constraint 
(just as in the continuum setting, see (jC4l) and (ICSp ). All three relate self-dual and anti-self-dual phase space 
parameters. 

In [2]| . we showed that provided we insist on a non-degeneracy condition,^ the constraint set consisting of 



It will turn out that these non-degeneracy conditions are crucial. In fact, the conditions do not only require 3d non-degeneracy 
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the matching constraints, the Gaufi constraints as well as the diagonal, cross and edge simplicity constraints 
{^<Lij,^(tij,(ti,&ij,&ijk,&ijki} allows one to construct an assignment of vectors yi{e) to the edges e of the 
triangulation. Here, the sub- index i indicates that j/i(e) is expressed in the frame of the tetrahedron {i}. The 
quadratic form of the constraints ([7]) allows for two possible sets. One assignment is such that the bi-vectors 
Xij can be written (modulo a sign) as the Hodge-duals of wedge products of vectors associated to two of the 
edges e, e' of the triangle {ij}: 

X,/=±*(y,(e)Ay,(e'))^. (8) 

Moreover, these edge-vectors are consistently transported between the frames of adjacent tetrahedra, that is 
yj°'{^) = Mji""^ yi'{e), where M""^ is the rotation matrix in the vector representation corresponding to the matrix 
]^AB bi- vector representation, see Appendix |B1 This latter property is ensured by the secondary or edge 
simplicity constraints and corresponds to the continuum condition that the spatial part of the spin connection 
is the Levi-Civita connection, i.e. it is determined by the spatial metric. 

In this so-called gravitational sector, the set of canonical variables {Xij,Mij} determines uniquely the as- 
signment ([S]) of vectors to edges modulo an orientation reversal for all vectors yi{e) i— > —yi{e). 

As we said, the constraints ([7]) allow for solutions in the so-called topological sector, which are of the form: 

X,/ = ±(2/,(e) A 2/,(e'))^. (9) 

These can be excluded using any one of two methods. On might replace some of the quadratic constraints with 
constraints cubic in the bi-vectors or one may add discrete conditions which fix the sign of certain quantities. 
This will be explained in more detail in Section [IV CI 

III. GAUGE- VARIANT ANALYSIS 

The constraints (O are scalar quantities and moreover, they are quadratic in the bi-vectors. Alternatively, 
the primary simplicity constraints can be recast in a format, which is linear in the bi-vectors. These conditions 
ensure that given the four bi-vectors Xij (varying {j}), one can reconstruct edge- vectors for the tetrahedron 
{«}. Thus, they have the advantage over the quadratic constraints (O that they kill the topological sector. For 
this reason, these constraints have been applied in the more recent construction of spin foam models [TT| - [l3| . 
This reformulation of the constraints requires an auxiliary vector, which if fixed to (1, 0, 0, 0), leads to what is 
known as the time-gauge. We shall adopt the time-gauge in our analysis and later detail how to move out of 
this gauge. 

If several tetrahedra share a given edge in the Sd-hypersurface, then we have also to guarantee that the vectors 
reconstructed in the different adjacent tetrahedra for that edge are identical, that is they are consistently parallel 
transported into each other. The secondary constraints, constructed below, see to this. They are the gauge- 
variant equivalent of the edge simplicity constraints and thus restrict the parallel transport matrices. As we 
shall derive later, they lead to constraints involving the bi-vectors of neighboring tetrahedra. 

A. Formulation of primary constraints 

Consider a tetrahedron Should one construct the bi-vectors as functions of the edge- vectors according 
to ([5]), then together they span only a 3d space, as a tetrahedron is a 3d object. Hence, for a geometrical 



but also that the normals of neighboring tetrahedra are non-parallel. In the particular case, where the 3d-hypersurface A is the 
boundary triangulation of a 4-simplex, these conditions ensure 4d non-degeneracy. 
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configuration, tliere exists a four-vector rii for every tetrahedron {i} such that: 

n,^ i*X,jr^ « 0, (10) 

for all four bi-vectors Xij in this tetrahedron. From now on, we shall refer to these gauge-variant conditions 
(especially in time-gauge, see below) as primary simplicity constraints. Note that this condition excludes the 
topological sector of solutions © (for non-degenerate tetrahedra). To see this, one should remember that the 
Hodge-dual -kX of a simple bi- vector X spans a 2d plane orthogonal to the plane spanned by X. Therefore, 
for a topological configuration, planes defined by the Hodge-duals of bi-vectors of the form ^ would share one 
direction, but three of these planes would span 4d space. 

On the other hand one can show ^3|] that bi-vectors satisfying the condition (flUl) also satisfy the diagonal 
and cross simplicity constraints. 

Were one to perform a gauge transformation on the reference frame of the tetrahedron z, the vector n.i would 
be rotated. We shall fix some of this freedom by choosing = et := (1,0,0,0) for every tetrahedron. This is 
known as time-gauge, and it yields a particularly tidy form of the simplicity constraints. We can formulate the 
primary constraints in the following manner: 

s,/ := {etT Xy "^ = - Xy■_o^ (11) 

where in this and subsequent formulae we shall reserve r, s, t to denote spatial indices taking values in {1, 2, 3}. 

Notably, the constraints become linear in the X variables but at the expense of introducing an explicit 
normal vector to each tetrahedron in the spatial hypersurface. A continuum version of this formalism has been 
developed in [3^ . 



B. Formulation of secondary constraints 

As stated earlier, now that one has reconstructed edge- vectors for each tetrahedron independently, one must 
ensure that these vectors agree for coincident edges. Thus, one must conceive of constraints on the connection 
degrees of freedom that map between the vector spaces associated to different tetrahedra. Let us consider 
two tetrahedra {i}, {j} with a shared face {ij}. Their reference frames have bases ^e" and ^e", respectively 
(a = 0,1,2,3). Moreover, we shall denote vectors by a; = Xi^ic"^. There exists an othogonal transformation 
between the two vector spaces: ^e" = Mji""- . Therefore, the components change like Xj°' = Mji"''^ Xi' . 

We are interested in the space of bi-vectors. We denote a generic bi-vector by Xi = Xi^ (^e A ^e)^, and 
(iC A ic)^ — ic"" A ic"" . The transformation between vector spaces induces a transformation between bi-vector 
spaces: X/ = M^/^ X,^ . 

A priori, we are not in a position to describe this rotation in the vector representation, since we do not 
start with the vector degrees of freedom. Rather, we start with the bi-vectors. The bi-vectors associated to 
the triangle shared by the tetrahedra {i} and {j} are denoted Xij and Xji. After the primary constraints are 
imposed, one is able to reconstruct the edge-vectors for each tetrahedron. We shall denote the edge-vectors 
associated to a particular edge in the triangle {ij} as yijk and Ujn, where the first index shall prescribe within 
which tetrahedron it is constructed. Put differently, we index an edge by {ijk} by a triple of tetrahedra sharing 
this edge. Moreover, by {jil} we usually refer to the same edge but expressed in the reference frame of the 
tetrahedron {j}. If an edge is shared by only three tetrahedra (as is the case for the boundary triangulation of 
a 4-simplex), one has k — I. 

To describe a geometric configuration, we must ensure that the three j/yfe per triangle are rotated into the 
three yjn . This will give us three constraints on the possible form of the matrix Mij . 

To be clear, a matrix in the vector representation is denoted by M"^, a matrix in the bi-vector representation 
is denoted by Af"^^, while once one moves to time-gauge, a matrix in the bi-vector representation is denoted 
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Our transformation may be separated nicely into three parts according to a number of criteria: 

X,,^ = -Afj/^ X,/ - - (3m,, ^Af,, (12) 

Note that *Xji — —Mji -k Xij if and only if Xji = —Mji Xij. Speaking in rough terms, any rotation can be 
decomposed into a part which preserves a plane pointwise, a part which rotates within that plane and a part 
which actively rotates that plane into another. 

Part 1 - ^Alji'. The first part of the SO (4) rotation preserves pointwise *Xij, that is: 

yr'' ^ CMJ^ y^r ^ y^'' if y^C*X,J. (13) 

Of the three parts, this is the one that is not completely determined by the bi- vector data and it is related 
to extrinsic curvature of the 3d-hypersurface. 

Part 2 - "^Mjii The next part rotates the plane of the triangle {ij} as seen in {i} into the plane as seen in {j}, 
that is, it rotates *Xij into *Xji: 

: Cm,, X,,f = -X,/. (14) 

Such a rotation Mji can be found as function of the bi- vectors Xij and Xji. 
Part 3 - ^Mjii Finally, we need the part which rotates the edge- vectors within the plane -ki^MjiXij): 

^M,, : (^Mj, y.juT ^ i^Mj, ^M,, y^.kT = (15) 

where {ijk} refers the edge shared by the triangles {ij} and {ik} and so on. The reason for the minus 
sign in the final equality will become clear when we define the edge- vectors. Once again, this rotation is 
can be determined as a function of the bi- vectors. 

All together, these three parts describe the most general transformation which ensures the geometricity of the 
3d-hypersurface triangulation. 

We shall now embark on making the above statements more explicit, beginning with: 
Part 1: We are interested in transformations which stabilise a particular bi- vector: 

X^ i^M X)^ = X"^ (16) 
There are two such transformations which satisfy this criterion: 

^M^"" ^exp{rnX-jf, and ^N^"" = exp {r^2 *X ■ , (17) 

where X := j^^. The J are the generators for S0(4), see Appendix [B] for precise definitions. Moreover, 
we are interested in the specific case that X is simple and that the transformations preserve, in a pointwise 
fashion, the plane defined by -kX. (This means that the vectors associated to the edges of the triangle are 
stabilised also.) The primary simplicity condition (in time-gauge) is given by (llip : 

(et)" e"''-^ X^ = 0. 

Equivalently, the condition states: X^'' = 0, where r,s E {1,2,3}. Thus, the transformations simplify to: 
^M^^ ^cxpimX^'-Jorf, and ^TV^^ = exp (r,2 1° V* Jst)^''- (18) 
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We can easily see that ^M"^^ preserves pointwise the plane defined by ★X, while it rotates the half- lines 

in the plane defined by X through an angle rji. On the other hand, ^N^^ preserves pointwise the plane 
defined by X, while it rotates the half-lines in the plane defined by X through an angle 772- Thus, the 
case we are interested in is that in which 772 = 0. 

At many points in our analysis, we split our degrees of freedom into self-dual and anti-self-dual variables: 

iM±^^=exp(,7±X±. Ji)^'', where r,± = (ryi ± r,2)&^. (19) 
Moving to time-gauge this yields: 

Vn±- := 2 ^M±0'-o« = exp ( ± x/2 ,± i± • j^)-, where ^ . ;= '^^^^0-*. ^''^ 

The ± sign arises since C'_|_0''0«C'gCOt _ gOrst 7 imposing simplicity and the preservation condition (772 = 0), 
we arrive at the conclusion: 

77+=77- = |, x+' = x.^ = X°^ (21) 
Thus, ultimately we have the condition that: 

= Cm-)-\ (22) 

To summarise, ^Mji^^ is a boost (in time-gauge). In other words, it is an S0(4) transformation effecting 
a rotation of the vector et, which is normal to the plane of the geometric triangle defined by -kXji. In 
time-gauge, this boost is mapped to two rotations mji± around the 3- vector x!j_ = x^_. The two rotation 
angles have the same absolute value but opposite signs. 

Part 2: Now, we are in the position to examine the part of the transformation which rotates the plane -kXij 
into -kXji. The transformation takes the form:* 

Xji± = -^Mji± Xij±, (24) 

where:^ 

n AB 



^Mji±^^ := exp 




— ' ' — IT ' 

Xji± X Xij-tll 



(26) 



We should note that ^M±^''os are not the only non-zero components of the matrix, but they are the only ones we need since we 
shall restrict to planes which lie in the 3d-hypersufa<;e perpendicular to et- 
* It is perhaps instructive to remember the transformation in 3d which rotates the vector a into the vector b: 

V = m^" a^, means that m = exp [ (sin~^ \b x a\)— — • e ) , (23) 

V |6xa| y 

where (er)st := trst and a := —. 
^ The expansion of an exponential in the adjoint representation is: 

exp(;U±B± • Ji)-*^ = cos(\/2/i±)P±'*^ -t- (1 - cos(V^Aj±))S±'^i±^ ± ^ sin(v^A'±)E±-^C±^^^ (25) 

v2 
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The transformation satisfying (|24p is not unique, as we can for instance multiply from the right any 
rotation leaving Xij invariant. However, this will not affect the arguments for the derivation of the 
secondary simplicity constraints we are going to make. 

In time-gauge ([26l) takes the form: 



2 {^Mj 



\OrOs 



exp 



(sm \xji± X Xij±\) 



Thus, here we can see that upon imposition of the simplicity constraint: 



(27) 



(28) 



In general proper rotations M, i.e. transformations leaving (1, 0, 0, 0) invariant, are mapped in time-gauge 
to two rotations m± with m+ = m_. 

Part 3: Ultimately, while this is enough to ensure that the triangles lie in the same plane, it is not enough to 
ensure that their edge- vectors are identified. For this, we need a rotation in the plane (^Af^ -kXij). To 
construct it explicitly, one needs to know the normalised edge-vectors. In this gauge-fixed context, the x 
can be viewed as the inward-pointing 3d-normals to the triangles within the hypersurface. Therefore, if 
the triangles {ij} and {ik} share an edge {ijk}, one can describe the edge- vector using the cross product 
between the two normals: 



yijk± 



\Xij± X Xik± I 



and 



yj^l± 



{Xji± X Xji±y 
\Xji± X Xji±\ 



In addition, we have already rotated the plane in which i/i^k is contained: 

yijk -> ( mji± yijk±) = ^ 



^ Xik± I 

Thus, the rotation in the plane which maps i^rriji yijk) into {—yju) is: 



\Xji± X Xik±\ 



exp 



/. 2 - ,\ {yjii± ^ mj^±yijk±) . 



\yjii± X "mjj± yijk±\ 



(29) 



(30) 



(31) 



The reason for the minus sign is due to the fact that edge-vectors are defined as the cross-product of 
the two normals in an ordering which necessitates its presence. Once again, we understand that once 
simplicity is imposed: 



1^ji,lk+ 



(33) 



The important point to notice about this part is that for a given triangle {ij} there are three ways to 
construct this rotation, which correspond to which pair of edge- vectors yijk, yju we choose. Obviously, 
these should all be equal since there is only one rotation which rotates the triangles into one another. We 
shall comment on this point shortly. 



In addition, it is rather easy to see what the gauge-non-fixed expression is: 



3 AB 



■■ exp 



where Yijk±' 



\\Xij± X 



1 \\Yjii± X ^Mji±Yijk±\\ \ *(Yjii± X ^Mji±Yijk±) j 



and so forth. 



(32) 
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With the exphcit form of the three rotations at hand, we can now determine the form of the secondary 
simphcity constraints. For every triangle {ij}, these must ensure that the rotation Mji consistently maps the 
three edge-vectors of the triangle {ij} from the reference frame {i} to the reference frame {j}. We constructed 
^Mji^Mji exactly such that is satisfies this requirement. Hence: 

{^Mj,''Mj,)-^Mj, (34) 

has to be a rotation that leaves the plane -kXij pointwise invariant, i.e. it must be of the form ^Mji. Therefore, 
it satisfies ([2^ in time-gauge. This leads to the secondary simplicity constraints: 

[{^mji^ik+ ^niji+y^ mji+ {^mji^ik+ ^niji+y^ mji-Y* = (35) 
which, since we are dealing with S0(3) matrices, can be encoded into the form:^^ 

where = function ({X}). 

The constraints (1361) involve quite complicated functions of the bi- vectors compared to the primary simplicity 
constraints (llip ). There are three constraints (with three components) per triangle {ij). These fix the rotation 
rriji- in terms of the rotation niji^, the rotation built out of a vectors Xij,Xji, and the rotation constructed 
from a pair of edge- vectors yijk, Vju- For the three possible pairs of edge- vectors, we obtain three conditions on 
the relation between rriji^ and rriji-. 

Note that the secondarysimplicity constraints resemble the reality condition for the complex Ashtekar vari- 
ables [U nil, (see also [23] for an alternative discretisation of these conditions). These reality conditions 
require that the real part of the spatial pull-back of the Ashtekar connection, i.e. the sum of the self-dual and 
anti-self-dual connections, is equal to (twice) the spin (Levi-Civita) connection. Here, {^mji^ik+'^mji+) is a 
discretization of the spatial spin (Levi-Civita) connection and for every triangle may be constructed out of the 
bi-vectors in three different ways - depending on the choice of edge in the triangle under consideration. To 
obtain an unambiguous discretisation of the spin connection, these three ways must coincide and this eventually 
leads to the so-called gluing conditions. 



Let us summarise the phase space variables and constraints in time-gauge. Our initial phase space parametri- 
sation in this gauge-fixed setting is given by the quadruple (xij±,Xji±,mij±,mji±) of data per triangle (ij). 
These variables feature the following Poisson algebra: 

{a;,,±°,m,,±''^} = ±2 

{x,,±°,mj,±'^''} = ±2 





7 


^ahc^cd 


7 


± 1 




7 


^ahc^dc 
^ '"'ji±^ 


7 


± 1 




7 


abc c 
c ■^ij± J 


7 


± 1 



{x,j±'',x^j±^} = ±2^^f''bc^c _ (37) 



This form of the constraints makes it explicit that there are at most three indpendent components per triangle and edge. But it 
allows for an ambiguity, namely the combination of rotation matrices appearing in the square bracket could also be a rotation 
with an angle tt. This has then to be excluded by hand and we shall implicitly understand this to be the case if one uses Il36p . 
(This problem does not appear if one uses the spin- 1/2 representation) . 
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Then, one imposes the constramt set: 



Sji Ik 



,rst I \st 

e (mij± mji±) 

3 



rst r/3 



Upon imposing the above constraints, one automatically satisfies the following condition: 

St 

Qji,lk±' 



^mjijk± {^mjij>k'±) ^ 



(38) 
(39) 
(40) 

(41) 
(42) 

(43) 



that is, the rotation in the plane is independent of the pair of edges chosen to construct the transformation 
^ji,ik- It turns out that one may replace some of the secondary simplicity constraints by the constraints (|43p . 



C. Relaxing time-gauge 



Time-gauge assumes that all the 4d-normals to the tetrahedra are rotated to et = (1,0,0,0). We can also 
formulate the simplicity constraints for the general case by incorporating the transformations needed to rotate 
the normals rii to et = (1,0,0,0). To start, we retain the splitting of degrees of freedom into self-dual and 
anti-self-dual variables x±'',m±^'' using the auxiliary vector e* = (1,0,0,0), see Appendix 1X1 

The primary simplicity constraints (|10p require that for every tetrahedron {i} there exists a (unit) vector n.i 
such that the Hodge-dual of the bi- vectors from this tetrahedron {i} are normal to this vector: n.i"' (■kXij)'^'^ — 0. 
We now consider an S0(4) rotation Ni"-'' rotating this vector rii to the standard one et — (1, 0, 0, 0): 

n- ^ TV,-" n," = et"". (44) 

We can split A^^ — (ni+,nj_) into self-dual and anti-self-dual components, so that the variables: 

Xij^ \ '^ij^ ^i— '^ij— (^^) 

and: 

my+™ = n,+^' my-+^* (^,+"')*", = n,.'"^ m,,_^* (^,-"')*" (46) 

are in time-gauge. Note that the rotation Ni is not uniquely defined, as one can multiply on the left by any 
rotation, which leaves et invariant. Hence, Ni should be understood as an element of SO(4)/SO(3). A rotation 
R leaving et invariant corresponds in the zb-splitting to two rotations r± satisfying r+ = r_. Thus, one can 
always choose Ni such that n+ = id, for example. The primary simplicity constraints can then be formulated 
in the following way (see for instance ^^): for every tetrahedron {j}, there exists an S0(3) rotation rn-^" such 
that: 

for all triangles {ij} in the tetrahedron {{]. Then, we define the variables Xij± and rhij± according to PSI 
For these variables, the secondary simplicity constraints (|^^ must hold, that is, we replace all the variables 
x± , m± there by their tilde counterparts x± , rh± : 

Sji,ik'' = e'^^* [{^mJ^^lk+ ^fhjt+y^ mji+ {^mjijk+ ^fhji+y^ f^iji^Y* . (48) 
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If we use n+ — id, we even have ^fhjijk+ '^fnji+ = '^n'>'ji.ik+ '^niji+, as here it is understood that the rotations 
'^'m'ji,ik+ and '^mji+ are fixed functions (|27[ I3ip of the variables a;+ — i+. The secondary simphcity constraints 
can then be expressed in a similar form to (j47p : 

Cmji^ik+ '^rriji+y^ mji+ {^mji^ik+ '^rriji+y^ w rii-mij- nj}. (49) 

The matching and closure constraints ((38l [39l |40)) can be expressed either in the x± , m± variables or in the 
x±,fh± variables, since these transform covariantly under rotations: 

^Cjji'' 2j^^.^r ^ grst (^f^^^^ mj,±Y^ = (n,± '^Z^j±Y (50) 

<± ^ Cji'' = Ej^»i±'' (nj± q±)''. 

Therefore, "'^c^i « if and only if ^tij± ~ and so on. 

The conditions (|47l49p involve the rotations nt- as auxiliary objects. To obtain constraints in the usual form, 
i.e. only containing phase space variables, one can expand the phase space to include the variables Ui- and 
their conjugate momenta. Then, one would have to add more constraints, which basically eliminate the Ui- 
and conjugate momenta as dynamical variables, so that in the end one obtains the same reduced phase space 
as before. Furthermore, the Gaufi constraints must be altered so that they also generate rotations for the ni_ 
variables. Work along these lines has been advanced in [s^, and plays an important role in the definition of 
the so-called projected spin networks [13 ■ A full canonical analysis is still open, however. 



IV. GAUGE-INVARIANT PHASE SPACE 



A. Gauge-invariant variables 



Here we shall discuss gauge-invariant variables in more detail. We shall be interested in investigating the 
relations among them arising purely from their definitions as well as the conditions imposed by the primary 
pT|) and secondary (|36| constraints. 

We shall define the gauge-invariant variables and discuss the relations among them using the time-gauge. The 
definitions for the gauge-invariant variables are also valid if we are not in time-gauge and just use the auxiliary 
vector et = (1,0,0,0) to define the S0(3) algebra (resp. group) valued variables x± (resp. m±) . The relations 
among the gauge-invariant variables and the restrictions imposed by the constraints are of course also valid off 
time-gauge. 

To start with, let us consider two gauge-invariant quantities that can be constructed from the data associated 
to one tetrahedron {i}: the squared areas aij± and the (outer) 3d-dihedral angles (j)ijk± between the normals 
Xij± and Xik± to the two triangles {ij} and {ik}, respectively: 



aij± 

cos (j)ijk± 



Xij±: ' Xij^ , 

Xij ' Xik 



^ Xij • Xij Xik ' Xik 



(51) 



There are four of the former and six of the latter for each chirality. In each ±-sector, the Gaufi constraints (|40 
give four relations among the six 3d-dihedral angles per tetrahedra: 



V«jj± + X! V"-ik± COS <pijk± = 



(52) 
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For each tetrahedron these Gaufi constraints can be solved for four of the 3d-dihedral angles, so that one is left 
with the four area variables and two (non-opposite) 3d-dihedral angles. This holds separately in the +- sector 
and the — sector. Hence, we obtain 2 x (4 + 2) = 12 independent gauge-invariant variables per tetrahedron. 
These are sufficient to determine the left-handed and right-handed geometry of each tetrahedron. 

From the 3d-angles (j), one can construct the (outer) 2d-dihedral angles a. This can be easily seen by 
considering directly the definition of the 2d-angles as products between the edge- vectors, for example yijk± and 

yiji±- 



cos atjki± 



(^Xij X Xif^ ) • i^X^j X Xil ) 



{Xij X Xik) ■ (Xij X Xik) \J {Xij X Xil ) ' {Xij X Xil ) 



COS (j)ikl - cos (pijk COS 7 

sin (jjijk sin (piji 



(53) 



We now turn to gauge-invariant data that require the consideration of more than one tetrahedron. Consider an 
edge {ijk} in the tetrahedron {i} and the same edge {jil} seen from the tetrahedron {j}.^^ We consider the 
angle between these two edges after parallel transporting {jil} to the reference system of {i}: 

[Xij X Xik) [Xji X Xjij . , 

cos0,k,ji± = ; / ' ' — . (54) 

Y [Xij X Xik ) ' [Xij X Xik ) Y [Xji X Xjl ) • [Xji X Xjl ) _j_ 

Hence, one can construct three angles, depending on the choice of edge, for every triangle. On a geometric 
configuration the three angles per triangle coincide and give the value of the (inner) 4d-dihedral angle between 
the two tetrahedra sharing this triangle (that is, the angle between the normal and reversed normal of the 
two neighboring tetrahedra) [2lj. Note that this interpretation does not necessarily hold for non-geometric 
configurations. We shall nevertheless refer to these angles as 4d-dihedral angles. 



Let us determine the dimension of the Gaufi-reduced phase space for a generic triangulation. We started with 
2x6 variables per triangle, hence the dimension of the initial phase space is 12 x Nt, where Nt is the number 
of triangles. In a 3d triangulation without boundary, every triangle is shared by two tetrahedra and every 
tetrahedron has four triangles. Thus, we have twice as many triangles as tetrahedra. This means that the 
number of Gaufi constraints totals 6 x ^ x A't. As these are first class constraints, we subtract twice their 
number from the dimension of the initial phase space. Therefore, the dimension of the reduced phase space 
is 6 x TVt. Remember that there are two area variables per triangle and four independent 3d-dihedral angles 
per tetrahedron, which amounts to four variables per triangle. ^'^ This leaves two further variables per triangle, 
in order to fully parameterise the Gaufi-reduced phase space. We must choose from the 4d-angles 6 defined in 
([5^ . of which there are three per triangle and chiral sector. 



This suggests that these six 6'-angles are not independent variables. Indeed, as we shall explain soon, they are 
related by terms involving the 2d-angles a. 



To analyse this interdependence, let us consider two tetrahedra {i} and {j}, as illustrated in Fig. [TJ We choose 
a gauge such that Xij± = (0, 0, 1) and Xji± = (0, 0, —1). Due to the conditions Xij± = —mij± Xji±, the rotation 



For the 3d triangulation defined by the boundary of a 4-simplex we have that I = k. 

We implicitly imposed the 'matching' constraints I I38I39I1 . They are second class constraints, but have trivial effect on the 
Poisson structure, and ultimately, all one needs to do is implicitly remember that if one writes down mij± and Xji±, that they 
are equal to {mji±)~^ and —mji±Xij±, respectively. 
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matrices mij± must now leave the vector — (0,0,1) invariant, i.e. they are S0(2) rotations. Thus, we can 
write mij± — R{S,±) where R{^) is a (positively-oriented) rotation around with angle ^. The angles f± are 
still not gauge-invariant quantities,^"* since there is a residual gauge freedom: one may still rotate around Cz 
in the frames of both tetrahedra. We can fix this residual gauge freedom by demanding, for instance, that the 
edge {ijk} in the tetrahedron {i} and the corresponding edge {jil} in the tetrahedron {j} are labelled by the 
normalised vectors yijk± = (0, 1, 0) and yjii± = (0,-1, 0), respectively. Here, we used again the definitions (p9)) : 

Xij X X'l}^ 

yijk± , , . , 

\J \Xij X Xik) ■ [Xij X Xik) _^ 

In this gauge, we obtain ■mij± — R(tt — 9ikji±). Note that the definition (IMl) for Oikji± does not fix the sign of 
the angle. This sign is determined by the equation mij± = R{tt — 0ikji±)- 





FIG. 1: Two tetrahedra sharing a triangle and the associated 2d dihedral angles and edge vectors. 



All the edge-vectors associated to the triangle {ij}, in the frames {i} and {j}, are in the plane spanned by 
Bx = (1,0,0) and Cy — (0, 1,0). Hence, one can write: 

yijk'± = R{aijkk'±) yijk±, and yjii'± = R{ajiw±) yjii±, (55) 

where again we augment the definitions (|53p of the a-angles so as to fix their signs. 

The angles Oikju Oik'ji' are attached to the triangle {ij}. They are computed using the pairs of edges {ijfc}, 
{jil} and {jik'}, {jil'}, respectively. For the relation between these angles we compute 

cos6ik',ji'± = yijk'± ■ 'mij± yjiii±, 

= [{R{aijkk') Vtjk) ■ R{t^ - 0ik.ji) {R{a]iU') y]ii)]± , (56) 
= [ilijk R{{Tr - Oikji - aijkk' + otjiiv) yjii]j^ . 



These angles are also used in [Tsl. l23l [2gt . 
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Let us reiterate that Oijkk' and ajiw correspond to the angle between the same two edges of the triangle {ij}. 
The difference is that aijkk' is computed from the data of the tetrahedron {i} whereas ajiw is computed from 
the data of the tetrahedron {j}. Let us define the sign of the angle 6ik'.ji'± in the same way as for 9ikji±. We 
finally obtain: 

Oik',jl'± — &ik,jl± + Ctijkk'± — 0:jiii'±. (57) 

Put into words, the ^-angles for a given triangle are defined using one edge but in two reference frames. Mean- 
while, any 2d-dihedral angle is defined using two edges but in the same reference frame. Equation ([Ff|) above 
states that the difference in (^-angles defined on two edges equals the difference in 2d-dihedral angles between 
those two edges as computed in the two reference frames. This explains why on geometrical configurations 
the three 0-angles per sector and triangle coincide: on these configurations the so-called gluing conditions hold 
[21I [27} . requiring that the 2d-dihedral angles at a triangle {ij} computed from the tetrahedron {i} should 
coincide with those computed from {j}. For our example, the gluing conditions impose that aijkk' ± = otjiiii±, 
which implies that 6ik'.ji'± = Oikji±. 



B. Restrictions imposed by the primary and secondary simplicity constraints 

Let us examine the implications of the primary (I4ip and secondary (|42p simplicity constraints for the scalar 
variables introduced in the previous section. 

The primary simplicity constraints (|4ip impose that Xij^ = Xij- for all triangles {ij}. We therefore conclude 
that any quantity, which can be computed without using the holonomies mij±, has the same value in both the 
self-dual and anti-self-dual sectors, that is: 

CLij-\- dij—., 

= (58) 
aijk+ = ctijk-. 

The secondary simplicity constraints (I42p are given by: 

[{^mijM+ "^rrii-j+y^ mij+ {^mij^ki+ ^mij+)"^ id- (59) 

In the gauge that we introduced in Section llV Al the matrices ^mij^ki+ and '^rriij^ are equal to the identity 
matrix. Consequently, we obtain the condition: 

mjj+ = {mij^y^ . (60) 

Since mij± = R{7r — 0ik,ji±) in this gauge, we can conclude (taking into account the 27r-periodicity of the angles): 

&ik,jl+ = —Qik,jl-- (61) 

Although obtained in a specific gauge, these equations hold generally as only gauge-invariant quantities are 
involved. In particular, we have the conditions (j6ip not only for the pair of corresponding edges {ijk}, {jil} but 
also for the other two pairs {ijfc'}, {jil'} and {ijk"}, {jil"}. Hence, we obtain three constraints per triangle: 

Qik,jl+ = ~(^ik,jl~ 

dik',jl'+ = —dik'Jl'- (62) 
(^ik",jl"+ = ^0ik" -. 
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Now, we are in a position to derive gluing conditions explicitly from the primary (j4ip and secondary (|42l) 
simplicity constraints. In the previous section, we derived that the 0-angles for the same triangle are related by 
differences in the 2d-dihedral angles a ([57)) : 

dik',ji'± — dikji± + aijkk'± — ctjiii'±. (63) 

If we substitute the conditions 9+ = — 0_ ([5^ and q+ = a_ ([55]) into equation we obtain the conditions: 

aijkk'± = ajiw±- (64) 

In other words, the 2d-dihedral angles belonging to a triangle, as computed from the two adjacent tetrahedra, 
have to coincide. These are the gluing conditions introduced and discussed in [2l|, [l^, HH . Here, we derived 
these gluing conditions explicitly from the primary and secondary simplicity constraints. 

The discussion above also allows one to derive an alternative version of the gluing conditions to those given 
in ([64)) . which uses the angles 6 instead of the 2d-angles a. Namely, along with the constraints (|64p . we also 
ensure the coincidence of the 4d-dihedral angles as computed from the three different edges of the same triangle: 

dik,jl± — &ik',jl'± — dik"Jl"± ■ (65) 

A similar relation can be derived for the 4d-dihedral angles of a (fiat) 4-simplex, where these 4d-dihedral angles 
are computed from the 3d-dihedral angles in different ways, see (27j, |3l|. Here, we generalised this relation to a 
phase space context, where the 4d-dihedral angles encode the extrinsic curvature. Hence, it is a free momentum 
variable and is therefore not determined by the 3d-dihedral angles. 

Note that whereas the primary simplicity constraints solely constitute conditions between the self-dual and 
the anti-self-dual sectors, the secondary simplicity constraints imply more. Specifically, they not only administer 
conditions between the sectors, but they also contain relations among quantities of one and the same sector. 



C. Relation to the diagonal, cross and edge simplicity constraints 

In this section, we shall reconsider the diagonal, cross and edge simplicity constraints ([7]): 

— Xijj^ ■ Xij^ — Xij^ ■ Xij^, 

— Xij+ ■ Xik+ — Xij^ ■ Xik-, (66) 
= Xik+ ■ (Afij+Xj;+) — Xik- ■ {Xij^Xji^), 

and show that these constraints follow from the primary and secondary simplicity constraints (|41|42p . 

To give ourselves a solid basis, let us restate some foundational facts. A geometrical configuration is one 
in which the bi-vectors are of the form X = *{y{e) A j/(e')), where y{e) denotes an edge-vector, and the 
holonomies ensure that the edge- vectors are assigned consistently across the tetrahedra. If the matching, Gaufi, 
diagonal, cross and edge simplicity constraints hold, one can prove 21] that on non-degenerate configurations, 
one has two possible ways to write the bi-vectors X^. They can be consistently written either as wedge 
products of vectors associated to the edges X^ = ±(j/(e) A y{e'))^ or as the Hodge-dual of such wedge products 
X^ = ± ★ (2/(e) A y{e'))^. The first case corresponds to the topological configurations, the second to the 
gravitational or geometrical configurations. The non-degeneracy conditions discussed in [21| do not only mean 
3d-non-degeneracy, i.e. the non-vanishing of 3d-volumes, but it also incorporate conditions pertaining to the 
embedding of the 3d-triangulation into the 4d- manifold. Namely, it prescibes that the normals of two neighboring 
tetrahedra must be non-parallel. (As an aside, remember that the normal to a tetrahedron is well-defined if 
the diagonal and cross simplicity constraints, as well as the 3d-non-degeneracy condition, hold.) Indeed in 
section|V]we will construct non-geometric configurations satisfying the Gaufi, diagonal, cross and edge simplicity 
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constraints, however featuring parallel 4d normals to the tctrahedra. This shows that the condition on non- 
parallel normals is necessary. 

Thus, if one shows that the primary and secondary simplicity constraints (|41I42I) imply the diagonal, cross 
and edge simplicity constraints then one has shown that these conditions are sufficient to ensure the 

geometricity of the triangulation for non-degenerate configurations. Furthermore, these constraints are also 
necessary, since they are satisfied by all geometric configurations. In fact, a nice feature of the primary and 
secondary simplicity constraints is that topological configurations are automatically excluded, as we shall see 
below. In contrast, since the simplicity constraints (j66p are quadratic in the bi-vectors, they can not distinguish 
between the topological and gravitational sectors, which are Hodge-dual to each other. In that case, to exclude 
the topological sector, one must add some conditions that are cubic in the bi-vectors (or replace some quadratic 
constraints by cubic ones, see also the discussion in [2TI|). 

Importantly, as we shall comment in Section |Vl the non-degeneracy conditions are essential if one uses the 
diagonal, cross and edge simplicity constraints (j66p . In contradistinction, the secondary simplicity constraints 
(|42p do not to rely on the non-degeneracy conditions and therefore constitute stronger conditions on the sector 
of 4d-degenerate configurations. 

This discussion about the equivalence of constraints (at least on the non-degenerate gravitational sector) will 
allow us to freely switch between the different constraint sets and to select the most convenient one in order to 
compute the Dirac brackets. This computation will be done in future work (3^ . To this end, we shall determine 
a set of independent constraints in the next section, which will also enable us to find the dimension of the 
reduced phase space. 

In Section IIVBI we derived conditions on the gauge-invariant variables from the primary and secondary 
simplicity constraints. We commence by showing that the diagonal, cross and edge simplicity constraints follow 
from these conditions. As shown in the previous section, the primary simplicity constraints imply (j58p : 

a,y4. = a,y _ and (j}tjk+ ^ 4>ijk- ■ (67) 

The diagonal and cross simplicity constraints in (j66|) are equivalent to:^^ 

aij+ = aij^ and cos (t>ijk+ ^ con (t>ijk- , (68) 

which clearly follow from (p7|) . The equations ([57)) also imply that the orientation in the -I— sector and — sector 
agree. To have a similar property with the quadratic constraints, one must add an extra constraint. For 
instance, one might require: 

C^^^e^^X,/X,,^X,,^ = 0, (69) 

which in time-gauge reduces to: 

Xij+ ■ {x.,k+ X Xii+) = Xij^ ■ [xik- X Xii^) (70) 
This states that the oriented volume squared as computed in the two sectors coincides. 

The secondary simplicity constraints imply (j6ip : 0y.fe/+ = —Oij.ki-- For a moment, consider the gauge-invariant 
quantity: 

Xik ' {j^ij Xjl ) 



cos Pik,jl± 



\J Xik ' Xik \JXjl ■ Xjl 



(71) 



^'^ We must assume non-vanishing areas. Moreover, for the rest of the discussion we shall assume 3d-non-degeneracy, i.e. the 
non-vanishing of areas, 3d-dihedral angles and 3d-volumes. 
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This is the 3d-dihedral angle between the triangles {ik} and {jl} (where the tetrahedra {i} and {j} share the 
triangle {ij})- Furthermore, there exists a relationship among various dihedral angles: 



cos Oik,jl± 



cos pikji - cos (pijk cos (pju 



sm (pijk sm (pjii 

Thus, if the 9ij^ki+ = -dij,ki- ^^^nd (j)ijk+ = 4>ijk- hold, then: 

cos pikjl+ = cos pikjl-. 

But the edge simplicity constraint, the last constraint in (|66p . appears in time-gauge as: 

Xik+ ■ {mij+Xji+) = Xik- ■ {rriij^Xji), 



(72) 

(73) 
(74) 



which is equivalent to ([75)) . 

As before, if one passes in the opposite direction, the edge simplicity constraint only implies that cos 9ikji+ = 
cos Oikji-- Therefore, one must add a further condition to the edge simplicity constraints, in order to fix the 
(relative) sign between the 4d-dihedral angles. Alternatively, one may consider a cubic version of the edge 
simplicity constraints [U [l^ : 

(75) 
(76) 



V A , ,r AB r<BCD Y C Y D n 

which in time-gauge reads: 

Xik+ ■ rnij+{xji+ X Xji+) + x^k- ■ mij^{xji^ x Xji^) = 0. 
This constraint implies Oij,ki+ = —Oij,ki~- 



D. Worked example: the boundary of a 4-simplex 

Our next task is to select an independent subset of simplicity constraints. This will facilitate a computation of the 
Dirac brackets in future work 34] and determines the dimension of the phase space corresponding to geometrical 
configurations. Here, we shall consider this question on the phase space reduced by the Gaufi constraints. Such 
a reduction simplifies hugely the calculation. Furthermore, we shall also restrict to the phase space associated 
to the 3d-boundary of a 4-simplex. More complicated triangulations will be the subject of further work. 

The boundary of a 4-simplex has 5 tetrahedra connected by 10 triangles. As explained in Section flV Al the 
Gaufi-reduced phase space is 60 dimensional and may be parametrised using gauge-invariant variables. The 60 
independent gauge-invariant variables needed to parameterise it are provided by the 20 area variables aij±, 20 
3d-dihedral angles (j)ijk± (two 3d-angles, associated to non-opposite edges, per tetrahedron and chiral sector) 
and 20 4d-dihedral angles 9ik,jk± (one 4d-angle per triangle and chiral sector). 

Now we move onto the constraint set. We shall show that the maximal independent subset of diagonal, cross 
and edge simplicity constraints contains 40 elements. We shall denote this set C. Given the parameterisation 
above, one can see that an independent subset of diagonal and cross simplicity constraints is given by taking 
all 10 diagonal simplicity constraints, and two cross simplicity constraints (acting on the angles associated to 
the same two non-opposite edges) per tetrahedron. These constitute the first 20 elements of C. 



Remember that I = k for the boundary of the 4-simplex. Also, note that one can explicitly reconstruct the other 3d-dihedral and 
4d-dihedral angles using I I521 1531 157| |. 
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We are left only to deal with the edge simplicity constraints. A priori, there are thirty edge simplicity 
constraints {cos 0ikjk+ = cos Oik, jk- or Oikjk+ = —Oikjk-), three per triangle. We shall argue that only 
twenty are independent. As illustrated in Section lIV A[ we can choose one constraint per triangle in (j62l) to 
explicitly relate quantities in the +-sector with corresponding quantities in the — sector. These make up the 
next 10 elements of C. We showed that the remaining edge simplicity constraints implicitly impose the gluing 
conditions ((64)) . which in this context take the form: 

0:ijkk'± = Ctjikk'±- (77) 

In words, the 2d-angles for a triangle {ij} coincide if computed from the two neighbouring tetrahedra {«} and 
{j}. Written in this form, there are 60 such gluing constraints, of which 10 are independent. We shall show that 
on the gauge-invariant phase space, with the 30 constraints already in C at our disposal, one has 50 relations 
among the gluing constraints ([77]) . 

For a start, note that the 2d-dihedral angles may be written in terms of the 3d-dihedral angles using the 
identity (|53l) . Using the diagonal and cross simplicity constraints in C, one has = a_. So if the gluing 
constraints hold for the +- sector, they will automatically hold for the — sector. We are left now with the 30 
gluing constraints on the +- sector, three per triangle. But remember on the gauge-invariant phase space, one 
has imposed the Gaufi constraint. This means that the tetrahedron is closed and furthermore it implies that 
the three 2d-angles for one triangle some up to tt. In this way, we obtain 10 more relations among the gluing 
constraints. 

The final 10 relations are quite involved since they are more non-local. But they can be understood to arise 
in the following manner. One may rewrite the gluing constraints as the conditions that the length of any given 
edge, computed in the three adjacent tetrahedra, must be unique. From these three constraints per edge only 
two are independent, that is we have another set of 10 relations. These are a manifestation of the final relations 
among the gluing constraints. Furthermore, these considerations can be confirmed by computing the Jacobian 
of the matrix of derivatives of the constraints with respect to the gauge-invariant variables [l^] . 

Thus, we have 10 independent gluing constraints, and we may pick them in a symmetric manner, say such 
that they simultaneously one per triangle and one per edge. These are the last 10 constraints to go into the set 
C. 

There is another way to formulate these last 10 gluing constraints. The 30 2d-angles a+ are determined, 
through the relations ([52)) and (|53)) . by the 10 squared areas a+ and the set of 10 3d-angles (/>+. We have 10 
independent gluing conditions between the 2d-angles. In light of the preamble, these descend to 10 independent 
relations between the areas and 3d-angles. For a geometric 4-simplex, the lengths of the 10 edges are free 
parameters. An alternative set of free parameters is provided by the areas. "'^^ Hence, the gluing constraints fix"'^^ 
the values of the 10 3d-dihedral angles 0+ as functions of the 10 (squared) areas: 

4>ijk+ = fijk{aim+)- (78) 

Unfortunately, an explicit expression for these functions fijk is not known. This problem is equivalent to 
expressing the ten lengths of a 4-simplex as a function of the ten areas, which however requires the solution of 
a polynomial of high degree with general coefficients. 

In summary, we end up with the set C of 40 independent constraints. C contains the 10 diagonal simplicity 
constraints and 10 cross simplicity constraints involving ±~pairs of two (non-opposite) 3d-angles per tetrahedron. 
These 20 constraints can be seen as primary constraints, ensuring the geometricity of each tetrahedron. Then, 
C contains another 20 constraints, which in a certain sense are conjugate to the first 20 constraints. We can 



This is the case because the matrix of derivatives of the areas with respect to the lengths is invertible at generic points, 
modulo discrete ambiguities 
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interpret these constraints as secondary as they ensure the consistent gluing between the tetrahedra. The 20 
constraints are formed from a set of 10 edge simphcity in the form co&9ikjk+ — cos6ikjk- and (10 further edge 
simphcity in the guise of) 10 gluing constraints: 4>ijk+ = fijk{o,im+)- 

These 40 constraints reduce the initial 60 dimensional (already Gaufi-reduced) phase space to a 20 dimensional 
one. This phase space can be parametrised by the 10 areas (or lengths) of the simplex as well as 10 4d-dihedral 
angles associated to the triangles of this simplex. The areas and the 4d-dihedral angles will be conjugated to 
each other [Hill]. 



V. DEGENERATE CONFIGURATIONS 



The diagonal, cross and edge simplicity constraints in the form ([7]) require that certain pairs of bi- vectors 
(including bi-vectors transported via a holonomy) span only a three-dimensional space. In particular, for two 
simple bi-vectors Xi"^ and satisfying: 

^AB ^ 0^ (7g) 

one can conclude that there exist three vectors u,v,w such that: 

X/ = (wAw)-^, X2^^{uAw)^. (80) 

For a proof, see for instance {l^ . 

But the simplicity constraints in the form ([7]) may fail to ensure geometricity on degenerate configurations. 



This was also noted in [21[ , where it was pointed out that the edge simplicity constraints also allow non-geometric 
configurations if the 4d-normals to two adjacent tetrahedra are parallel. For instance, consider a configuration 
on the boundary of a 4-simplex, where all the bi-vectors -kX are orthogonal to the unit vector (1,0,0,0). In 
other words, we are in time-gauge. Additionally, let us say that the SO (4) matrices M"^^ are, in fact, elements 
of its S0(3) subgroup that leaves (1, 0, 0, 0) invariant. Then, all the bi-vectors entering the simplicity conditions 
are orthogonal to (1,0,0,0) and so any pair of these bi-vectors spans only a three-dimensional space and the 
diagonal, cross and edge simplicity constraints are satisfied. Importantly, however, this is not a geometric 
configuration. 

In fact, note that this subspace of configurations corresponds to an S0(3) BF phase space. A related 
space of configurations was found in a saddle point analysis for the EPR spin foam model in [ssj . For the 
example taken, namely, the boundary of a 4-simplex, there are 30 independent gauge-invariant quantities 
in this subspace of the phase space. This is the dimension of the GauB-reduced S0(3) BF phase space. 
Moreover, it is 10 dimensions larger than the subspace containing the geometric configurations. For the subspace 
containing degenerate configurations, the dihedral angles 6, as defined in cannot be interpreted as 4d- 
dihedral angles. This interpretation is only valid for geometrical configurations, while the condition that these 
degenerate configurations satisfy is that the information they encode is confined to a 3d space. 

This shows that if one uses the form ([T]) of the constraints, one must also require 4d-non-degeneracy. In 
particular, the normals of neighbouring tetrahedra should not be parallel. This non-degeneracy requirement 
appears also in every spin foam model [131, see in particular the discussion in [s^. Essentially, one considers 
there the configuration space associated to a 4-simplex: the 10 so(4)-valued bi-vectors associated to the triangles 
of the simplex, which satisfy the diagonal and cross simplicity constraints as well as the Gaufi constraints. The 
connection is assumed to be flat. For non-degenerate configurations, one can show that these requirements allow 
the reconstruction of the 4-simplex geometry. Put differently, one can consistently assign vectors to the edges 
of the 4-simplex [s^]. This geometry is characterised by 10 gauge-invariant quantities. We may choose the 10 
lengths or the 10 areas of the 4-simplex. For degenerate configurations, however, one obtains a 15-dimensional 
space of gauge-invariant quantities, see also 
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Once again, the edge simplicity constraints in the form (O fail for 4d-degenerate configurations. We wish to 
point out that this is not the case for the gauge- variant form of the secondary simplicity constraints psp , which 
therefore constitute a stronger requirement concerning 4d-degenerate configurations. 

Also, recasting some of the edge simplicity constraints as gluing conditions aijkk' = ajifcfc' restricts the 
4d-degenerate configurations. The gluing conditions do not completely exclude degenerate configurations, but 
they do however reduce the dimensionality of the degenerate sector. In effect, one its left with degenerate 
geometries, while non-geometric configurations are removed. Therefore, they could also be applicable to spin 
foams where they might suppress the degenerate sector by reducing its dimensionality. 

Whereas for one 4-simplex, the reduction eliminates 5 dimensions and still leaves the areas (encoded in 
the spins associated to the triangles) as free variables, the situation is much more complicated for a larger 
triangulation - both for the 3d hypersurfaces and the 4d triangulations. Then, for geometric configurations, the 
areas are not any more free variables. Actually, there are (non-local) constraints among the areas to ensure that 
they follow from a consistent assignment of edge-lengths in the triangulation (whose number is typically smaller 
then the number of triangles) |38|| . The formalism of area-angle Regge calculus [273 allows one to express these 
constraints in a local manner by keeping the 3d-dihedral angles as free variables. Then, the gluing constraints 
([M)) together with the GauB constraints lead to restrictions on the areas of the triangulation. 

To summarise, the condition of non-degeneracy plays a crucial role in the reduction from topological BF 
theory to the gravitational sector. For spin foam models, it would therefore be important to study whether the 
dynamics encoded in the amplitudes and measure passes interchangable among non-degenerate and degenerate 
sectors. Were the answer affirmative, then one would have to devise methods to impose non-degeneracy condi- 
tions on the measure and amplitudes. Possibilities for this include the gluing constraints (j64p and the secondary 
constraints in the form ([35]). 



VI. DISCUSSION 



Spin foam models attempt to provide a path integral for gravity by starting from a partition function for 
topological BF theory and implementing simplicity constraints within the amplitudes and measure. The cor- 
rect imposition of these constraints is therefore central to the entire approach. In current spin foam models, 
one typically deals only with the primary simplicity constraints, but it has been pointed out that both 
primary and secondary simplicity constraints are necessary in order to produce the correct expectation values. 
Additionally, there should be a relationship between the boundary Hilbert space defined by spin foam models 
p^ . [20| and the one which follows from a canonical quantisation, where one would have to deal with both types 
of constraint. 

Therefore, we studied different formulations of the primary and secondary simplicity conditions in this work. 
We have seen that these different formulations interact differently with the degenerate sector that plays an 
important role in the quantum dynamics described by spin foam models (35| . 

In particular, we derived gauge-variant simplicity constraints. In this form, the primary constraints (|4ip 
are linear in the bi- vectors. It is in this guise that they recently led to proposals for new spin foam models 
[lll,[l^. The secondary constraints ([32]) involve the holonomies and, through an explicit construction of the 
spatial Levi-Civita connection, the bivectors in a complicated way. These constraints resemble discrete versions 
of the so-called reality conditions for (complex) Ashtekar variables, which ensure that the sum of self-dual and 
anti-self-dual connections is equal to (twice) the spatial Levi-Civita connection. 

A well-defined spatial Levi-Civita connection requires one to be on a proper geometric configuration, i.e. the 
gluing conditions (|64p must hold. For this reason, the gluing constraints actually follow from the secondary 
simplicity constraints (or, equivalently, from the edge simplicity constraints in the non-degenerate sector). The 
gluing conditions can be alternatively expressed as constraints (|65p ensuring that the extrinsic curvature, i.e. 
the 4d-dihedral angles are well-defined. This can be easily understood by realising that the self-dual or anti- 
self-dual connections are a linear combination of the spatial Levi-Civita connection and the extrinsic curvature. 
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Hence, an ambiguity in the definition of tlie Levi-Civita connection ieads to an ambiguity in the definition of 
the 4d-dihedrai angles and vice versa. 

From just the primary simplicity constraints one might expect that the effect of the simplicity constraints 
is to reduce the S0(4) phase space to an S0(3) phase space. This is however not the case, since the gluing 
conditions (which are contained within the secondary simplicity constraints) reduce the SO (3) phase space even 
further. If one uses the diagonal, cross and edge simplicity constraints (O, it is essential that one requires 
non-degeneracy during this part of the reduction,. On the other hand, if one uses the secondary simplicity 
constraints in the form (|^^ and in the form of the gluing constraints ([M]) . one automatically removes the 
degenerate sector. A further analysis of this issue would serve to benefit the assembly of methods to exclude 
non-degenerate configurations from spin foam models. 

In this work, we formulated different constraint sets and elucidated the relations between these various sets. 
We discussed the dimensionality of the reduced phase space for the case of the boundary of a 4-simplex. For 
more complicated triangulations, the issue is much more involved [2]| as one has to deal with non-local area 
constraints. Another issue is the derivation of the symplectic structure, i.e. the Dirac brackets for the reduced 
phase space. This will be the subject of further work [sj], where we shall also comment on the role of the 
Immirzi parameter in the different stages of reduction. 



Appendix A: Conventions 



The discretisation of 4-dimensional SO (4) BF-theory relies heavily on the isomorphism between the space of 
bi- vectors (anti-symmetric second order tensors) and the algebra so (4). In our conventions the commutation 
relations among the generators of the algebra take the form: 

= J^. (Al) 

The indices are A = aa where a,a £ {0, 1, 2, 3}. The structure constants are: 

qABC gArs gBst gCtr 

^ab ^dbC ^ab ^abC ^ab ^abC ^ab gabC 

where the Kronecker delta and the summation convention on this space are given by: 

gAB gab^ab gabgab^ 

j^AB^BC Ij^Abb^bbC so that ^AB^BC^^AC^ 

We can use the Hodge ★ operator to split the algebra into self-dual and anti-self-dual sub-algebrae: 



(A2) 



(A3) 



J±^:=F±^^J^, where P±^^ := i(<5-4^ ± e^^). (A4) 

It is straightforward to check that these projectors are orthonormal, that is P^^ P^'~^ = Sss'Pf'^ , furthermore 
P^^^ + p AB _ gAB j^^y then proceed to generate all manner of projected quantities: 

- P±^^' X^' , M±^ = P±^^' M^' , (A5) 

(j_^ABC _ Pj_^^' (jA'BC _ p^BB' qAB' C _ p_^CC' qABC' _ p_^AA' p^BB' p_^CC' qA' B' C' ^^g^ 

The following identity for the 50(4) structure constants: 

^^ABC C^A'B'C ^ 2\P±^^' P±^^' - P±^^' P±^'^). (A7) 
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is in close analogy to the relation e^^'^e"''^''' = {5'"'' 6^''' - S"^'S"'^) for the 50(3) structure constants. 

We can explicitly map from the self-dual and anti-sclf-dual sectors to corresponding S0(3) quantities in the 
vector representation. To this end, we use an auxiliary unit vector, which in this work is Cf = (1,0,0,0). We 
then define: 



m±''^ = 2 et« et* P±«''-^ P±*«<^ M^<^. 



where r, s take values in {1, 2, 3}. 



Appendix B: Various representations of SO(4) 

Let us directly construct the relationship between various representations of S0(4) for future reference: 

Vector: The generators of the vector representation are: {J"^)^ = —S^^, and thus a generic transformation is 
given by: 

M"'' = exp(S • J)"^ (Bl) 

Bi-vector: One can then generate the bi-vector representation by taking the tensor product of two vector 
representations and antisymmetrising: 

M^^ :=2M[«[6M»l5] =exp(-B- J)^^, where {J^f^ = C'^^^. (B2) 

Quaternion: The vector representation is isomorphic to the representation carried by the quaternions H. One 
maps a 4-vector x = {x° , x^ , x"^ , x^) to a pseudo-unitary matrix x = x^ao +ix^ar.^^ Then, the action of 
S0(4) on this space is: 

-1 u ^± ■= exp(-|6±-CT), 

X Q+xq , where (B4) 
^ b±- := B±°\ ^ ' 

Adjoint: Furthermore, the bi-vector representation is isomorphic to the adjoint representation. One maps a 
bi-vector to an element of so (4): X = X^J^. Then, SO (4) acts on its algebra by conjugation: 

G := exp(-iS-J), 

X ^GXG-^ =G+X.(G+)-^ +G_X_(G)-\ where , (B5) 

G± := exp(-iB±. J±). 
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Appendix C: Continuum canonical formulation of Plebanski theory 

We describe some properties of an initial canonical analysis of the continuum theory. In [l^, where a 
comprehensive analysis of the Plebanski action is undertaken, they use a slightly different form of the simplicity 
constraints. It has the same class of non-degenerate solutions. Here, upon performing a Legendre transform, 
one finds that the phase space is parameterised by: {wi^,Ili^ := \ eoijk CP-y ^jk)^), where V^^^ := 
i((5"^^ + ie"^^). The parameters wq"^ and (p^^ are Lagrange multipliers imposing the primary constraints: 

D,n,A ~ 0, 

(5^C ^BC _ 1 ^AB ^CD^ x^C n,)^ « 0, ^^^^ 

where Di is the covariant derivative with respect to the spatial components of the spin connection. The first 
constraint is the Gaufi constraint and generates S0(4) transformations, while the final constraint is the canonical 
form of the simplicity constraints. One should then check the preservation of the constraint hypersurface under 
time evolution. This yields secondary constraints. We shall be interested in the secondary constraint arising 
from the preservations of the simplicity constraint above: 

IS^C gBC _ 1 ^AB ^CD^ Xo..C J^^^^.^(^_^2)Di^, ^^^^H ^^BS ^^R (^_^ n,)^] , (C5) 

where C^^^ are the structure constants of so (4). We note that this constraint is a functional of both the 
bi-vector and the connection. We shall see a reflection of this in the discrete theory. 
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